Spin noise in polariton lasers 



m 
o 

(N 

P-h 

(N 



i 



i 

C 

o 
o 



> 
m 

(N 
O 
m 

> 

X 



M.M. Glazov, 1,2 M.A. Semina, 1 E.Ya. Sherman, 3 ' 4 and A.V. Kavokin 2 ' 5 

1 Ioffe Physical- Technical Institute of the RAS, 194021 St. -Petersburg, Russia 
2 Spin and Optics Laboratory, St. -Petersburg State University, 198504 St. -Petersburg, Russia 
3 Department of Physical Chemistry, The University of the Basque Country, 48080 Bilbao, Spain 
J> IKERBASQUE Basque Foundation for Science, Bilbao, 48011 Bizkaia, Spain 
5 School of Physics and Astronomy, University of Southampton, 
Highfield, Southampton S017 IB J, United Kingdom 

We develop a theory of spin fluctuations of exciton-polaritons in a quantum microcavity under the 
non-resonant unpolarized pumping. It is shown that the spin noise is sensitive to the scattering rates 
in the system, occupation of the ground state, statistics of polaritons and interactions. The spin 
noise spectrum drastically narrows in the polariton lasing regime due to formation of a polariton 
condensate, while its shape becomes non-Lorentzian owing to interaction-induced spin decoherence. 
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Introduction. Quantum microcavity is a system where a 
semiconductor quantum well is placed between the Bragg 
mirrors making it possible to achieve the strong coupling 
between the light and matter. In this system the energy 
is coherently transferred back and forth between the pho- 
ton trapped in the microcavity and the exciton, the ele- 
mentary excitation of the semiconductor. First observed 
by Weisbuch et al. [Tj, the strong coupling results in the 
formation of mixed light-matter quasiparticles, exciton- 
polaritons, extensively studied since then [21 [3J. The 
exciton-polaritons combine an extremely small effective 
mass, inherited from the photon, and strong interactions 
between themselves and with the environment due to the 
cxcitonic fraction in this quasiparticle, allowing one to 
study plenty of fascinating phenomena including polari- 
ton lasing [3], which consists in accumulation of a huge 
number of polaritons in a single quantum state (non- 
equilibrium condensate) and spontaneous emission of co- 
herent light by this state. 

The exciton-polaritons are characterized by the spin 
projections s z = ±1 onto the structure growth axis z cor- 
responding to the right or left circular polarization of the 
photon and to the same spin component of the exciton. 
Superpositions of s z = ±1 polariton states give rise to 
the linear or elliptical polarization of exciton-polaritons. 
Polariton lasers represent a model bosonic system where 
the polarization of light emitted by a microcavity directly 
corresponds to the spin state of exciton-polaritons, allow- 
ing one to study the polariton spin dynamics by optical 
methods [5 . A wealth of prominent spin-related effects 
were proposed and realized in microcavities, including 
self-induced Larmor precession [51 [7] , linear polarization 
inversion [FJ, optical spin Hall effect [9HTT]. spin Meiss- 
ner effect [T2"HT4"j . spin multistability [13 US], etc., see 
Refs. [17] for reviews. 

In polariton lasers, the spontaneous symmetry breaking 
results in the appearance of a stochastic vector polariza- 
tion including circular polarization directly proportional 
to the spin of a polariton condensate [TSJ [TH] . The Stokes 
vector (pseudospin) of a condensate may be pinned to one 
of the crystal axes [20], due to the structure anisotropy. 
In this case the time-averaged polarization of emission of 



polariton lasers is defined by pinning, while it momentary 
value may fluctuate. The spin fluctuations give rise to the 
"spin noise" which may be studied experimentally by op- 
tical means. The spin fluctuations, being inherent to any 
system at equilibrium or not, were first observed in 1981 
in Na vapor [21]. Nowadays they became a research field 
in semiconductor spintronics, since measurements of their 
properties can provide an important information on spin 
dynamics, hardly accessible otherwise [22H2I]- 

Here we study theoretically the spin noise spectra of 
polariton lasers both below and beyond the laser threshold 
and demonstrate that it is extremely sensitive to the mean 
occupation number of the condensate and to the statistics 
of polaritons. Our approach can be easily extended to a 
rich variety of bosonic systems where spin polarization 
may be formed including cold atoms, indirect excitons, 
triplet superconductors and 3 He superfluids. 

Model. The dynamics of the exciton-polariton spin dou- 
blet can be described with the pseudospin approach where 
the density matrix of polariton state with the wavevector 
k can be written as p k — N k I + Sk ■ tr. Here I is the 
unit 2x2 matrix, <x is a pseudovector composed of Pauli 
matrices, N k is the spin-average occupancy of the state 
fe, and Sk is the pseudospin of polariton in this state. In 
what follows we focus on the dynamics of the ground state, 
corresponding to k = and treat all other states in the 
system as a reservoir. Correspondingly, we omit k sub- 
script in the notations. The pseudospin components S a , 
where a = {x,y, z} is the Cartesian index, correspond to 
the polarization of emission: S z /N gives its circular po- 
larization degree, while S x /N and S y /N give the linear 
polarization degree in the axis frames (xy) and (x'y') ro- 
tated by 45° with respect to each other. [2] The dynamics 
of the occupation number and the polariton pseudospin is 
governed by the set of kinetic equations [28 
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where tq is the lifetime of polaritons in the ground state, Q 
is the effective magnetic field arising from (i) anisotropy 
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of the system and (ii) from polariton-polariton interac- 
tions, its explicit form will be given below. The scalar, 
Q n {N, S}, and pseudovector, Q S {S, N}, collision inte- 
grals describe arrival and departure of the particles into 
the ground state. The general form of the collision in- 
tegrals is presented in Ref. [35] for the polariton-phonon 
interaction and in Ref. for the polariton-polariton 
interaction. Here we adopt their simplest form 
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W out N -W m {l + N), (3) 

(w out - w in + r s - x )[s - s (si)]. (4) 



Here W out and W m are out- and in-scattering rates, re- 
lated with the presence of the reservoir, as schematically 
illustrated in Fig. [I] In particular, W m is proportional 
to the occupation of the reservoir and it is determined 
by the pumping rate. The term proportional to 1 + N 
reflects the bosonic nature of exciton-polaritons. It de- 
scribes the stimulated scattering processes. In the colli- 
sion integral for the polariton spin, Q S {S, TV}, t s is the 
spin relaxation time and So(f2) is the steady-state spin 
induced by the effective field SI. If polariton exchange 
with reservoir is efficient, one can introduce the effective 
temperature T of the polariton system. The steady-state 
spin is -(N)Sl/\Sl\ for T < fl and for T > SI. In what 
follows (unless otherwise is specified) , we consider the case 
of high temperatures, where (S) = 0, and the occupancy 
of the polariton ground state given by the balance of in- 
and out- scattering processes: 



W out 
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where the condition W m < t^ 1 + W out should hold [30] , 
The fluctuations of the condensate occupation number 
5N(t) ee N(t)-{N), and pseudospin 5S{t) = S(t)-{S) = 
S(t) are described by the correlation functions, namely: 
£(t) ee (5N{t')SN(t' + t)) and C a0 (t) ee (S a (f)S p (t' + tj), 
where the angular brackets denote the averaging over the 
time t' for a given t. According to the general theory of 
fluctuations [3"TH3"3"] the correlation functions of the fluc- 
tuations as functions of t or t' obey the same set of kinetic 
equations as fluctuating quantities [34 . The solution of 
Eqs. and ^ in the absence of effective magnetic fields 
and interactions (SI ee 0) , results in the exponential time- 
decay of correlations and isotropic spin fluctuations: 

£(t) = /C(0)e-l*l/ T % C aP (t) = 6 af! C aa {a)e-W T °<*, (6) 

where 8 a p is the Kronecker <5-symbol, single time correla- 
tors (mean square fluctuations) JC(0) — ((SN) 2 ), C aa (0) = 
(S 2 )o will be found below, while the particle-spin corre- 
lations (6N(t')S a (t' + t)) vanish and will be disregarded. 
We introduced the correlation times r c , t CiS according to 
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Equations (|6| and Q clearly show that the particle num- 
ber and spin fluctuations of exciton-polaritons decay ex- 
ponentially as a function of time difference \t\ and the 
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FIG. 1: (a) Illustrative scheme of the pumping. The reservoir 
and the ground state are shown as well as in- and out- scatter- 
ing processes, (b) Temporal dependence of spin fluctuations 
calculated for the (TV) = 10 (green/solid), 100 (blue/dotted) 
1000 (black/dash-dotted). Other parameters are: to = 25 ps, 
W out — 0, t s = 10 ns. Dashed curve is calculated for t s — > oo 
and (N) — 1000. (c) Spin noise power spectra calculated for 
the same parameters, (d) Spin noise power spectra calculated 
with allowance for the anisotropic splitting Sl a = 2 ns -1 . Other 
parameters are the same as in panels (b), (c). 



correlation time of the particle number fluctuations is r c , 
while the spin fluctuations are correlated on a shorter 
time scale r CtS . The spin fluctuation spectra defined as 
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C a p(t)e lu)t dt have a Lorentzian form: 
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with the half-width half maximum t~} determined by 
the inverse spin correlation time. It follows from Eq. ^ 
that the correlation time of fluctuations is strongly en- 
hanced and the spin noise spectrum is strongly narrowed, 
if (N) 3> 1, i.e. where the ground state is macroscopi- 
cally occupied. In particular, if r s — > oo and W out —> 
which corresponds to the negligible spin-flip in the ground 
state and negligible depletion of the condensate due to 
the scattering of polaritons back to the reservoir, we have 
for the correlation time of particle and spin fluctuations 
T c — r o(l + (A0)- The life-time of exciton-polaritons in the 
state-of-the-art structures varies from ~ I ps to ~ 100 ps 
with the corresponding maximum (N) being between 10 3 
and 10 5 , yielding r c in the range of 1 ns ... 10 /is. Corre- 
spondingly, the spin noise frequencies range from MHz to 
GHz for the macroscopically occupied ground state. The 
typical spin noise spectra and temporal dependence of the 
spin fluctuations are presented in Fig. [ljb),(c) where the 
drastic effect of the ground state occupation is clearly vis- 
ible. The several orders of magnitude enhancement of 
the spin correlation time and corresponding narrowing of 
the noise spectrum is a purely bosonic effect, inherent to 
exciton-polaritons. Indeed, any fluctuation (in the par- 
ticle number or in the spin) is supported by the bosonic 
stimulation due to 1 + N factor in the in-scattering term 
of the collision integral in Eq. Q . 
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The mean square of the particle and pseudospin fluc- 
tuations can be found using the approach of Ref. |35j 
(see also Ref. [3]) where the statistics of the exciton- 
polariton ground state was found by means of the mas- 
ter equation approach. Assuming that interactions and 
effective magnetic field are absent, the system becomes 
spin-isotropic and can be described by the independent 
occupations of two orthogonal spin states, N-f and N±, 
described by the same distribution functions P(A^). Us- 
ing these distribution functions one can express the mean 
square fluctuations, that is t — correlators ((5N) 2 ) = 
Eiv tlJ v, P(W t )P(^)[(tfr + - (iV)] 2 and (S 2 a ) Q = 

En„n, P(A T )P(A i )[(A T - A ; )/2] 2 , as 

(Sl)o = ((SN) 2 ) = \(N)[1 + (g^ 1)(A)], (9) 

where (S^)o corresponds to isotropic fluctuations, and 
g^ is the second order coherence of a single state. In 
particular, g^ = 2 corresponds to the thermal statis- 
tics |36j . where particle and spin square fluctuations are 
oc (A)(l + (A)) and grow quadratically with the ground 
state occupation. This situation is realized for equilibrium 
Bose gas [37] or if the ground state feedback on reservoir 
is negligible. By contrast, g( 2 ~> — 1 corresponds to the 
coherent statistics, in which case the fluctuations are sup- 
pressed, being cx (A). In the limit of low occupancy of 
the ground state, ((8N) 2 ) 1 / 2 = \J (A) in agreement with 
the theory of classical gas. Moreover, the full statistics 
of spin fluctuations can be determined, as a convolution 
of distribution functions P(N - S a ) and P(N + S a ). We 
present here the result in the (A) ^> 1 limit: 

Pcoh(Sa) = MA))" 1 / 2 exp (~S 2 J(N)), g^ = l, (i ) 
p th (S a ) = (N)- 1 e W (-2\S a \/(N)), g^=2. (11) 

For pulsed excitation, the stochastic polarization of po- 
lariton condensates was studied in Ref. |19j by numeri- 
cal integration of Langevin equations for the condensate 
spinor wavef unctions. The approach of Ref. [JJ5] can be ex- 
tended to a continuous wave pumping making it possible, 
in principle, to determine the full statistics of the system. 
Here we resort to a simpler analytical approach based on 
kinetic equations treating g^ as an external parameter. 

Role of effective magnetic fields. We begin the discus- 
sion of the effective magnetic fields with the case of an 
anisotropic system where the polariton doublet is split 
into the pair of states, linearly polarized along x and y 
axes. In such a case the vector f2 = (f2 a ,0,0) deter- 
mines the anisotropic splitting 6 , 38 - 40] , and we assume 
that the effective temperature T of the system exceeds the 
anisotropic splitting to neglect the steady spin polariza- 
tion. The S z and Sy show precession, while the dynamics 
of S x remains purely dissipational. As a result, the fluc- 
tuations become anisotropic with nonzero spectrum com- 
ponents (c.f. Ref. [41] ) 
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where C 



C&\ = [C%>]*, (S 2 a ) is given 

by Eq. j9j, and r C:S by Eq. Q. Figure Wd) shows 
the spin noise spectra calculated with the allowance for 
the anisotropic splitting for different occupations of the 
ground state. It is seen that the single peak is trans- 
formed to the two peak structure. This is because the 
spin correlation time t c s depends strongly on the ground 
state occupation: for small pumping rates and small (A) 
the product Q a T c,s ^ 1 an d the splitting is not visi- 
ble, however, for larger pumping n a r C)S becomes com- 
parable or larger than 1, making the anisotropic split- 
ting Ml a resolvable in the spin noise, as demonstrated in 
Fig. |TJd) . If the temperature is so low that T <C fifl a , 
the polaritons are predominantly polarized along the ef- 
fective field direction x. In this case, (S x ) — —(A), 
(S y ) = (S z ) =0, while the mean square fluctuations take 
the form ((S x - (S x )) 2 ) = ((SN) 2 ), (S 2 y ) = (S 2 ) = (N)/2. 
Moreover, the non-trivial single time correlation appears: 
(SyS z ) — (S z S y )* = —i(N)/2. Spin noise spectrum in 
this case can be obtained in a similar way. Specifically, 
Eqs. ( fl2| ), ( [13] ) still hold, but with different numerators 
corresponding to mean square fluctuations found above. 

Now let us consider the effect of polariton-polariton in- 
teractions on the spin noise. As follows from multiple ex- 
perimental and theoretical works [HI H2T - TI41 [TBI |4"2"] these 
interactions are strongly spin-anisotropic: the exciton- 
polaritons with the same z pseudospin components, i.e. 
with the same circular polarizations, repel efficiently due 
to the exchange interaction of electrons/holes with the 
same spin, while the polaritons with opposite circular po- 
larizations can weakly attract each other [SI H2] , the latter 
effect is neglected here. Hence, the interactions create an 
effective fluctuating field f2 = (0, 0, fij), directed along the 
z axis. Here the magnitude Mli = a\S z (otx > 0) is related 
to the fluctuations of polaritons pseudospin z component. 
This effective field has two important consequences on po- 
lariton spin dynamics and spin noise: (i) it induces preces- 
sion of the pseudospin around the z axis, known as self- 
induced Larmor precession [S] [7J, and (ii) it suppresses 
fluctuations of the z pseudospin component, favoring lin- 
ear polarization of the macrooccupied state [TU H3J . 

It is instructive to start the analysis with the spin pre- 
cession effect assuming that the effective temperature is 
high enough, T » a\^J (S\)q. We also ignore in what fol- 
lows the anisotropic splitting, taking fl a — 0. Clearly, the 
effective field Hi results in the dephasing of the in-plane 
pseudospin components. To describe the effect quantita- 
tively we note, that according to Eq. ^ the correlation 
functions of the in-plane pseudospin components are gov- 
erned by the following equation: 

^# + ^±^(t)S,(t) = 0. (15) 
at t c s n 

Here the upper (lower) sign corresponds to a = x (a — y), 



^xx 



2t„ 
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correlations between z— and in-plane components are dis- 
regarded because the effective field i~2 does not couple 
them, and, under the above assumptions, S z can be con- 
sidered as an independent parameter whose fluctuations 
are given by Eqs. ^6|, ffify. The set of linear equations ( 15 1 
can be readily solved as: 



= e -\t\h 



exp 



.Oil 



1*1 



(16) 



with (. . .) z meaning the averaging over the fluctuations of 
S z . The treatment of the general case is beyond the scope 
of the present work, here we consider two limiting cases: 
the regimes of (a) fast and (b) slow fluctuations. 

If ai(S z )^ 2 T CiS /h <c 1, the change of the effective field 
takes place on the time scale, which is much smaller than 
the spin precession period in the interaction induced field 
fli. This case corresponds to the motional narrowing and 
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exp 
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r, 2 



ksi) T CtS \t\ 



(17) 



In this limit the spin fluctuations decay exponentially, re- 
sulting in the Lorentzian shape of the spin noise power 
spectrum. This spectrum is, however, anisotropic: the 
half-width half maximum of and C y ^ is larger than 
that of dfz because the interaction-induced field does 
not affect S z owing to the spin-anisotropy of polariton- 
polariton interactions. 

In the opposite limit, where the fluctuations of spin z 
component are slow enough, so that the in-plane pseudo- 
spin components make several turns during the correlation 
time of S z , i.e. if ai(S z )^ 2 T c , s /h ^ 1, the fluctuations 
of S z and, hence, the effective field Hi can be considered 
frozen. In this case, the spin dephasing takes place on 
the timescale of . The particular t— dependence of 
the correlator C aa (t) is determined by statistics of the 
condensate [44] . For (N) > lwe obtain for two important 
limiting cases of coherent and thermal statistics (t <§; t c s ): 



Caa (*) 
Caa. (^) 

Mo 



= exp(-r c 2 oh * 2 ), 



(J 



(2) _ 



1 



.(2) 



2. 



(18) 
(19) 



where the dephasing rates are: r 2 oh = a\(N)/Ah 2 and 
r 2 h = a\{N) 2 /-ih 2 = (N)r 2 oh . In this limit, the temporal 
dependence of the spin fluctuations is directly related to 
the ground state statistics: Gaussian fluctuations of S z 



described by p co h(S z ) in Eq. (10 1 result in the Gaussian 
decay of the in-plane spin components, while the sharply- 



peaked exponential distribution pth(S z ) in Eq. (11) for 
thermal statistics, results in the slow power-law decay of 
the fluctuations due to high probability of small values 
of S z , corresponding to low precession rates. As a result, 
the noise spectrum of the in-plane pseudospin components 
deviates strongly from Lorentzian. Making Fourier trans- 
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FIG. 2: (a) Temporal dependence of the in-plane pseudospin 
fluctuations calculated with allowance for interactions ac- 
cording to Eqs. (|TsJ) and (19} for (TV) = 1000, «i = 
10 -2 ns -1 [42] for coherent statistics (red/solid) and thermal 
statistics (blue/dotted), (b) Spin noise power spectra calcu- 
lated for the same parameters as in (a). 



form of Eqs. (18 1, (19) we obtain: 



^act 

Wo 



■ exp 



(-c 2 /4r 2 oh ), 9 ^ = l, (20) 



ih 



■exp(-H/T th ), g 



(2) 



(21) 



It is noteworthy, that the non-analytical dependence of 
the spin fluctuations for the thermal statistics at uj — > is 
related with the power-law temporal decay of the fluctu- 
ations. The typical temporal dependence of the in-plane 
pseudospin fluctuations and noise power spectra are plot- 
ted in Fig.[2](a) and (b), respectively. Figure clearly shows 
different qualitative behavior of the noise of the in-plane 
pseudospin components for different statistics of the po- 
laritons. We emphasize that the drastic difference of the 
decoherence times and noise spectral widths for coherent 
and thermal statistics is related to different dependences 
of (S z ) on the ground state occupancy, (N), see Eq. ^. 

Note, that if ol\ (S z ) is comparable with or larger than 
the effective temperature T of the system, the fluctua- 
tions of the pseudospin z component would be suppressed. 
This effect can be modeled by multiplying the distribution 



function of S z in Eqs. ( 10 1 and ( 11 ) by the Boltzmann fac- 
tor exp (—haiS 2 /T) for the probability of thermal fluctu- 
ations [T2] . As a result, in the limit of T — > we obtain for 
(S z ) = T/2hcti. In this case, the dephasing rate of the in- 
plane pseudospin components, which determines the spin 
noise spectral width, can be estimated as ~ a/ a\T jh. 

Conclusions. A theory of spin fluctuations in polari- 
ton lasers had been developed. We have demonstrated 
that the spin noise of exciton-polaritons strongly depends 
on the occupation numbers, statistics and interactions in 
the polariton condensate, thus providing a tool for the 
studies of a large variety of bosonic spin effects. Various 
regimes of spin noise have been identified. Experimental 
verification of these predictions can be done by Fourier 
spectroscopy of Kerr of Faraday rotation spectra. Exten- 
sion of this model to other spin-polarised bosonic systems 
is straightforward. 
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